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摘要 : 星 前 核 是 研究 恒星 形成 初始 条 件 和 恒星 形成 极 早期 阶段 的 重要 样本 。 密 度 和 温度 是 描述 
星 前 核 最 基本 的 两 个 物理 参数 ， 也 是 利用 非 局 域 热 动 平 衡 辐 射 转移 计算 方法 研究 分 子 激 发 条 件 
[0 计算 分 子 柱 密度 ， 进 行 天 体 化 学 模拟 计算 所 必需 的 物理 参数 。Barnard 68 (B68) 是 一 个 非常 
典型 的 星 前 核 ， 它 具有 类 似 球形 的 结构 ， 非 常 适合 用 来 检验 球 对 称 的 理论 模型 。 在 球 对 称 的 坐标 
系 中 ， 利 用 多 方 的 气体 静 力 学 模型 拟 合 了 B68 的 尘埃 消光 数据 ， 模 型 给 出 的 星 前 核 B68 的 质量 
为 0.8 Mo。 与 之 前 等 温 的 气体 静 力 学 模型 相 比 ， 该 模型 不 仅 能 拟 合 尘埃 的 消光 数据 ， 且 拟 合 结 
果 与 B68 中 气体 非 等 温 分 布 的 观测 结果 相符 ， 另 外 与 压强 的 观测 数据 也 符合 得 很 好 。 

关 8 词 : 分 子 云 核 ， 恒星 形成 ， 星 前 核 ， 流 体力 学 
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1 引 言 


分 子 云 中 的 致密 核 是 恒星 形成 的 场所 。 星 前 核 (prestellar core) 是 自 引力 束缚 的 致密 分 
子 云 核 。 星 前 核 可 以 在 自 引 力作 用 下 南 缩 形成 原 恒 星 ， 并 最 终 形 成 恒星 呈 。 星 前 核 代 表 了 整 
个 恒星 形成 过 程 中 最 初始 的 状态 和 极 早 期 的 阶段 ， 对 于 设 定 恒 星 形成 数值 模拟 和 理论 研究 
的 初始 条 件 提供 了 很 强 的 约束 。 因 此 ， 星 前 核 是 研究 恒星 形成 初始 条 件 和 极 早 期 阶段 的 
重要 样本 。 

星 前 核 的 物理 条 件 ， 如 密度 和 温度 ， 可 以 通过 尘埃 的 消光 、 人 尘埃 的 发 射 和 分 子 谱 线 的 观 
测 进行 测量 。 例 如 ， 通 过 人 尘埃 消光 的 观测 数据 可 以 得 到 星 前 核 的 柱 密度 分 布 ， 结 合 模型 模拟 
进而 获得 星 前 核 的 密度 结构 ， 而 通过 远 红 外 和 亚 毫米 波 的 多 波段 观测 数据 可 以 同时 获得 星 
前 核 的 密度 结构 和 人 尘埃 的 温度 结构 ”"。 星 前 核 的 气体 温度 需要 通过 分 子 谱 线 的 观测 获得 ， 
tn NH". 

利用 尘埃 和 分 子 谱 线 的 观测 数据 所 估算 出 的 密度 和 温度 结构 ， 是 描述 星 前 核 最 基本 的 
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物理 参量 。 在 非 局 域 热 动 平 衡 的 辐射 转移 计算 中 ， 密 度 和 温度 是 研究 分 子 谱 线 激 发 条 件 、 模 
拟 分 子 谱 线 轮廓 、 估 算 分 子 丰 度 进而 研究 分 子 耗 散 和 人 尘埃 冰 形 成 的 必要 输入 参数 ”“"。 另 


外 ， 在 进行 恒星 形成 理论 研究 >“”"、 数 值 模拟 ”和 天 体 化 学 模拟 计算 ”的 过 程 中 ， 密 度 和 温 


度 是 最 基本 的 物理 量 。 
Barnard 68 (B68) 


是 一 个 典型 


的 星 前 核 ， 具 有 近 


圆 形 结构 ( 见 图 加， 可 以 用 来 检验 恒星 


形成 相关 的 模型 "。 利 用 光学 和 近 红 外 波段 的 观测 数据 ，Avles 等 人 "精确 地 测量 了 星 前 核 


B68 中 的 侍 埃 对 背景 星 的 消光 ， 通 过 观测 数据 构 开 
( 见 图 四 b))。 人 尘埃 的 消光 可 以 通过 N(H+2H2)/Av==2x10?! cm~?-mag 


了 尘埃 消光 随 着 星 前 核 半 径 的 径 向 分 布 
-1 转换 成 氧 分 子 


的 柱 密 度 ， 进 而 可 以 研究 B68 的 密度 结构 。 进 一 步 研究 发 现 ，B68 的 密度 结构 可 以 很 好 地 


示 ， 尤 其 是 分 子 谱 线 的 观测 数据 表明 ，B68 的 气体 温度 并 非 等 温 , 
左右 ， 边界 大 概 6 ~ 7 K”。 因 此 ， 有 必要 利用 更 精确 的 模型 
的 密度 结构 和 温度 结构 。 本 文 的 第 2 BEATE 


用 等 温 的 静 力 学 模型 ( 即 等 温 的 莱恩 - 埃 姆 登 方 程 ) 来 描述 ”。 但 是 ， 后 来 更 多 的 观测 结果 显 


步 又 ;第 3 章 为 模拟 的 结果 和 讨论 ; 第 4 章 对 全 文 进行 总 结 与 展望 。 
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光 的 径 向 分 布 中。 


2 模型 与 拟 合 


2.1 ”尘埃 消光 数据 


RA(J2000) 
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YE: a) 数字 巡天 (Digitized Sky Survey) 的 光学 图 像 ， 


而 是 中 间 温 度 高 为 10 K 
新 研究 B68 这 一 典型 星 前 核 
多 方 气体 静 力 学 模型 ， 以 及 进行 模型 拟 合 的 基本 


红色 的 “十 ”号 表示 持 埃 消 光 最 大 的 位 置 ，b) 尘埃 消 


图 1 星 前 核 B68 的 光学 图 像 和 尘埃 消光 的 径 向 分 布 


这 里 所 使 用 的 星 前 核 B68 的 数据 来 自 参 考 文献 由。 利用 光学 和 近 红 外 波段 对 B68 区 域 


星光 的 观测 ，Avles 等 人 ”精确 地 测量 了 星 前 核 B68 中 的 尘埃 对 背景 星 的 消光 ， 得 到 了 人 尘埃 


消光 的 径 向 分 布 。 
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2.2 ”多 方 的 气体 静 力学 模型 

星 前 核 B68 可 以 看 作 一 个 球 对 称 的 理想 气体 球 ， 并 且 在 引力 和 热 压力 共同 作用 下 处 于 
流体 静 力学 状态 。 由 于 B68 的 气体 温度 分 布 为 内 部 高 边界 低 ”， 因 此 为 了 更 合理 地 描述 
B68 的 密度 结构 和 温度 结构 ， 我 们 采用 球 对 称 的 多 方 气体 静 力 学 模型 只 。 质 量 密度 为 : 


P = po” = [gH Nnu, 0” , (1) 


其 中 ，po 和 ny, 分别 是 星 前 核 中 心 的 质量 密度 和 数 密度 ，7nas 是 H 的 质量 ，He 是 权重 ， 
是 多 方 指数 。0 由 下 面 的 微分 方程 所 决定 ， 即 : 


1 d /dN on 
ere, ee © 


其 中 ，& 是 一 个 无 量 纲 的 半径 ， 与 星 前 核 半径 之 间 的 关系 为 : 


iy) 
z 


r= y (n+1) KoU" /(4rG)E ， (3) 


其 中 ，K 是 常数 ，G 是 引力 常数 。 星 前 核 中 心 的 质量 密度 应 该 是 有 限 的 ， 式 回 的 边界 条 
TEA: 


6(0)=1, 9(0)=0 . (4) 


是 


前 核 气体 的 压强 为 : 


P=Kp\t= , (5) 


根据 理想 气体 的 状态 方程 ， 利 用 式 (四 和 (6) 可 以 得 出 星 前 核 的 气体 温度 结构 为 : 


_ WmaP _ umy 
kg p kg 


EP, my 是 氧 原子 的 质量 ，ks 是 玻 尔 兹 曼 和 常数 ， 是 平均 分 子 权 重 。 设 星 前 核 的 中 心 温度 
NT, 由 式 (0) 可 以 得 出 : 


Ke 1/n (6) 


Toke Tok 
K=— I/n 一 : jn * ”) 
LM Po HMH (ugm, nH, ) 
星 前 核 的 质量 为 : 
M=-—4n | (n + 1)Tokg/(umy) a ia sg (8) 
| 4nG ae &=fmax | 


RP, Emax 为 最 大 的 无 量 纲 半径 值 。 通 过 这 个 多 方 的 气体 静 力 学 模型 可 以 研究 星 前 核 B68 
的 密度 结构 、 温 度 结构 、 压 强 结构 并 计算 其 质量 。 


[Some 


ChinaXiv 合 作 期 刊 


224 天 文学 进展 40 卷 


2.3 ”模型 拟 合 

利用 多 方 的 气体 静 力 学 模型 拟 合 B68 的 尘埃 消光 数据 (图 四 b))， 其 基本 步 又 如 下 所 述 。 
首先 ， 给 定 一 个 多 方 指 数值 ， 结 合式 团 ， 利 用 四 阶 的 龙 格 - 库 塔 数值 计算 方法 可 以 求 出 式 
(2) 的 一 个 数值 解 。 其 次 ， 给 定 星 前 核 的 一 个 中 心 数 密度 ny, 值 ， 可 以 得 出 星 前 核 的 一 个 数 
密度 分 布 为 ntsb"(6)。 第 三 ，B68 中 心 的 气体 温度 为 二 =10K””， 把 上 述 的 n, ma A To 
的 数值 带 入 式 (四 并 结合 式 ( 国 ， 可 以 得 出 云 核 半径 x 和 & 的 关系 。 第 四 ， 给 定 一 个 € 的 数 
值 ， 就 可 以 得 到 星 前 核 的 半径 。 另 外 ， 由 于 观测 到 的 尘埃 消光 数据 是 角 尺 度 的 函数 ， 所 以 需 
要 给 定 星 前 核 的 距离 D， 从 而 把 模型 给 出 的 线 尺度 转换 成 角 尺 度 。 最 后 ， 把 得 到 的 数 密度 随 
着 半径 变化 的 分 布 沿 着 视线 积分 得 到 柱 密度 ， 通 过 N(H+2H2)/Av =2x10?! em?-mag-1™ 
把 柱 密 度 转 换 成 尘埃 消光 ， 然 后 对 比 理论 上 的 消光 数值 与 观测 的 尘埃 消光 数据 。 综 上 所 述 ， 
在 整个 计算 过 程 中 ， 模 型 包括 n, na, 上 和 D 共 4 个 自由 参数 。 为 了 寻找 与 观测 数据 符合 最 
好 的 参数 值 ， 我 们 采用 蒙特 卡 洛 优化 算法 ， 通 过 优化 软件 包 MAGIXZ@ 实现 ” 。 


3 ”结果 与 讨论 


es 利用 蒙特 卡 洛 的 优化 算法 ， 得 到 与 观测 数据 相符 的 理论 模型 ( 见 图 团 。 图 四 给 出 了 等 温 
© 的 气体 静 力 学 模型 与 多 方 的 气体 静 力学 模型 拟 合 观测 数据 的 对 比 结果 。 相 较 于 等 温 模型 ， 该 
多 方 模型 给 出 了 更 大 的 半径 。 图 田 显 示 了 模型 中 4 个 自由 参数 的 后 验 函数 分 布 ， 这 些 一 维 的 
(© 后 验 函数 非常 接近 高 斯 分 布 。 分 布 函数 的 最 大 似 然 值 为 最 佳 拟 合 参数 值 ， 误 差 取 自 分 布 函数 
Q 的 16% 和 84% KA HA BT MEWS BBA AS Te n = 8.02111, Æ 

前 核 的 无 量 纲 半径 € = 2.46492, uty Ho 的 数 密 度 nu, = 3.1649 x 108 cm-3， 以 及 星 


前 核 的 距离 D = 74.5805! pes 


A,/mag 


10 100 
R/C") 
注 ， 带 有 误差 棒 的 圆 点 为 实测 的 尘埃 消光 数据 ， 黑 色 实 线 为 等 温 的 气体 静 力学 模型 只 ， 蓝 色 实 线 为 本 文 所 采用 的 


多 方 的 气体 静 力 学 模型 。 


图 2 ”观测 数据 与 理论 模型 的 对 比 


° https: //magix.astro.uni-koeln.de/ 
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n=8.0241.9 


E=2,4670-2 


-0.14 


ny /10°=3. 16735 


n, /10° 


D=74.58 0 


3.0 3.06 318 3.30 735 765 
n é n,, /10° D 


YE: 最 优 参 数值 取 自 最 大 似 然 值 ， 左 右 误差 分 别 为 16%% 和 84% 的 分 位 数 。 


图 3 “利用 蒙特 卡 洛 优化 算法 得 到 的 参数 的 一 维和 二 维 后 验 函 数 分 布 图 


归 一 化 


--- 压强 分 布 
10 100 
R/(") 


YE: 中 心 的 温度 和 压强 分 别 为 10 K 和 5.14x1071! Pa. 


图 4 归 一 化 的 模型 温度 和 压强 随 着 云 核 半 径 的 变化 分 布 


226 


把 拟 合 的 最 人 


天 文学 进展 


ChinaXiv 合 作 期 刊 


40 卷 


参数 值 带 入 式 ( 团 和 (@), HAAR (M) 可 以 得 出 星 前 核 B68 的 温度 结 


构 和 压强 结构 。 图 团 显示 了 归 一 化 的 温度 和 压强 的 径 向 分 布 ， 其 中 心 的 温度 和 压强 分 别 为 
10 K 和 5.14x10-1 Pa， 边 界 的 温度 和 压强 分 别 为 6.7K 和 1.3x107! Pa。 相 较 于 之 前 的 等 


温 模 型 ”"， 该 模型 给 4 


考虑 星际 辐射 场 


泡 (Loop 1 superbubble) 的 附近 。 
圈 1 超 泡 的 压强 为 1.2x10-12 Pa 左右， 该 模型 得 
近 观 测 值 。 最 后 ， 通 过 式 (四)， 可 以 得 出 
量 最 终 转化 成 恒星 质量 ”， 在 


M57”, 


与 之 前 Li 等 人 RERA, RIEA 
百 星 结构 的 基本 方 
F 也 比较 苛刻 ， 必 须 存在 一 个 径 向 的 线 习 


知 ， 是 研究 ， 
模型 ， 满 足 的 条 伯 


JRS 


常 严格 的 限制 。 


4 总 结 


sO) =H 


在 球 对 称 坐 标 系 中 ， 我 们 采用 多 方 气体 静 力 学 模型 对 B68 HA 
在 观测 数据 与 模型 的 拟 合 过 程 中 ， 利 月 
数 ， 进 而 得 到 了 最 佳 的 拟 合 参 
气体 静 力 学 模型 不 仅 能 拟 合 观 涡 
模型 给 出 的 星 前 核 B68 的 质量 为 0.8 Me， 假 设 气 体质 量 转化 成 恒星 质量 


与 展望 


前 核 B68 


的 边界 温度 值 更 符合 6 ~ 7 K 的 观测 结果 一 。 与 气体 的 温度 测量 结 
果 相 比 ， 边 界 的 尘埃 温度 测量 值 更 高 ， 这 主要 是 由 星际 辐射 场 造成 8 ”于 。 
影响 ， 因 此 模型 无 法 解释 尘埃 的 温度 测量 结果 。 另 外 ，B68 位 于 圈 1 超 
伦琴 X 射线 天 文 台 (ROSAT) H X 射线 观测 结果 可 知 ， 
的 边界 压强 为 1.3x10-1 Pa， 非 常 接 
的 质量 为 0.8 Mo。 假设 30% 的 气体 质 
自 引 力 的 作用 下 有 望 形成 一 个 0.24 Mo 的 恒星 ， 其 光谱 型 为 


自 引 力作 用 下 有 望 

阿 塔 - 
到 银 道 面 上 距离 较 远 的 星团 
型 可 能 会 有 一 定 的 局 限 拆 
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The Polytropic Static Gas Model of Prestellar Core 
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Abstract: Prestellar core is important to improve our understanding about the initial condi- 
tions and the earliest stage of star formation. Density and temperature are the basic physical 
parameters to characterize the prestellar core, and also are two required input physical pa- 
rameters to study excitation conditions and estimate the column densities of molecular lines 
under non local thermodynamic equilibrium, and run the astrochemical models. Barnard 68 
(B68) is a typical prestellar core and its structure is quite similar with a sphere, which makes 
it a proper source to test star formation models. Under the spherical coordinate system, 
the extinction data of B68 is fitted using polytropic static gas model, which gives a mass of 
0.8 Mo for B68. Compared to the previous isothermal static gas model, the model adopted 
here not only can fit the dust extinction data very well, but also agrees with the observational 


non-uniform gas temperature, and furthermore can meet the observational pressure. 


Key words: molecular cloud core; star formation; prestellar core; hydrodynamics 


